Dynamic disorder in receptor-ligand forced dissociation experiments 
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Recently experiments showed that some biological noncovalent bonds increase their lifetimes when 
they are stretched by an external force, and their lifetimes will decrease when the force increases 
further. Several specific quantitative models have been proposed to explain the intriguing transitions 
from the "catch-bond" to the "slip-bond". Different from the previous efforts, in this work we 
propose that the dynamic disorder of the force-dependent dissociation rate can account for the 
counterintuitive behaviors of the bonds. A Gaussian stochastic rate model is used to quantitatively 
describe the transitions observed recently in the single bond P-selctin glycoprotein ligand 1(PSGL- 
1)— P-selectin force rupture experiment [Marshall, et al., (2003) Nature 423, f90-f93]. Our model 
agrees well to the experimental data. We conclude that the catch bonds could arise from the 
stronger positive correlation between the height of the intrinsic energy barrier and the distance from 
the bound state to the barrier; classical pathway scenario or a priori catch bond assumption is not 
essential. 

PACS numbers: 87.15.Aa, 82.37.Rs, 87.15. By, 82.20.Uv 



Since Bell firstly proposed the famous force induced 
dissociation rate [l|, 

fc off = fc° ff exp[^], (1) 

where fc° ff = fco exp(— f3AG$) is the intrinsic rate con- 
stant in the absence of force, AG* is the height of the 
intrinsic energy barrier, x* is projection of the distance 
from the bound state to the energy barrier along the 
applied external force /, and = k^T with fee the 
Boltzmann's constant and T absolute temperature, the 
expression has been demonstrated experimentally 0, Q 
and widely employed in various forced dissociation ex- 
periments. Later, at least four other models have been 
put forward to explain and understand newcome forced 
dissociation experiments [3- 0- E3| • Inparticular, Dembo 
proposed a Hookean spring model 0, @ to describe force 
responses of receptor-ligand bonds. In addition to pre- 
dicting that the dissociation rates of the bonds increase 
exponentially with the square of the force, the most im- 
portant contribution of the model may be the finding of 
a "catch bond" which is defined as increasing its lifetime 
when the bond is stretched by the force. Correspond- 
ingly, a bond described by Bell expression is defined "slip 
bond" for its lifetime decreases when the force is apply- 
ing. 

Until recently, the catch bond predicted mathemati- 
cally was demonstrated in some biological adhesive bonds 
which include the lectin-like bacterial adhesion protein 
FimH 0, P-selctinglycoprotein ligand 1(PSGL-1)-P- 
or L-selectins 0, 1 1 1| complex. However, these experi- 



ments also observed that the catch bonds always transit 
into slip bonds when the stretching force increases be- 
yond a certain values, i.e., their lifetimes are shortened 
again. The counterintuitive catch-to-slip transition has 
attracted considerable attention from experimenters and 
theorists. Several kinetic models have been proposed to 
explain the intriguing observations in qualitative and 
quantitative approaches 0, 0, Q . We know that the 
interface between the ligand and receptor in the adhe- 
sive complex has been reported to be broad and shallow, 
such as the crystal structure of PSGL-1— P-selectin com- 
plex revealed [l^. In addition, as one type of noncova- 
lent bonds, the interactions between the molecules are 
weaker. Therefore it is plausible that the height and po- 
sition of the energy barrier of the complex fluctuate with 
time due to thermal motion of the whole macromolecu- 
lar structure. Dissociation reactions with fluctuating en- 
ergy barriers have been studied in terms of rate processes 
with dynamic disorder , which was first proposed and 
theoretically investigated by Agmon and Hopfield [l7| . 
Hence it is of interest to determine whether the fluctu- 
ation of the height and position of the energy barrier 
induces the catch-to-slip transition. On the other hand, 
we also note that in Bell's initial work and in the other 
models developed later, the intrinsic rate constant fco and 
the distance x* were determined and time-independent. 
It is possible to derive novel results from relaxation of 
this restriction. Stimulated by the two considerations, in 
the present work we propose a stochastic Gaussian rate 
model to quantitatively describe the catch-to-slip bond 
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transitions. In addition to well predicting the experi- 
mental data, our model provides a new possible physical 
origin of the catch bonds: they are likely to be induced 
by stronger positive correlation between the fluctuating 
height of the energy barrier AG* and the distance ar. 

Consider a simple molecular dissociation process under 
a constant force /, 



Binding state (B) — > Unbinding state (U), 



(2) 



where the time-dependent forced dissociation rate kf(t) 
is a stochastic variable. If the survival probability P(t) 
of the state B is assumed to satisfy the first order decay 
rate equation then its formal solution is given by 



P(t) = (exp(- j o k f (T)dr)\ . (3) 
Cumulant expansion of the above equation |l8| leads to, 



AG* and x^ are used, we name it Gaussian stochastic rate 
model (GSRM). 

The average dissociation rate Eq. JHJ is so simple that 
we can immediately distinguish four different physical 
situations according to the definitions of the parameters: 
(i) if all K x , K g , and K gx vanish, then average rate is 
just the classical Bell expression Eq. (Q; (ii) if both K x 
and K gx vanish or in the absence of the fluctuation of 
the distance x*(t), (kf(t)) still keeps the Bell formula 
except that the intrinsic dissociation rate changes into 
ko exp(—PAG^+(3 2 K g /2), i.e., the fluctuation of the bar- 
rier height speeds up the dissociation process [17$. (iii) 
If both Kg and K gx vanish or in the absence of the fluc- 
tuation of AG* , then we have 



(k f (t)) = fc exp^AG*) exp 



P4f 



(3 2 K X 



(7) 



P(t) w exp 



— / d,T(kf(T)) + higher orders terms 



(4) 



where the higher order terms mean the higher correla- 
tion functions of the rate k(t), e.g., the double integral 
of the second order correlation function (kf(ti)kf(t2)) — 

(*/(*i))<*/(*2)> etc. 

According to standard Arrhenius form, we rewrite 
the Bell expression as (kf(t)) = fc Q (exp[— /3(AG*(i) — 
fx$(t)))). Now the characteristic of the Bell expression is 
determined by the two stochastic processes AG^i), the 
energy barrier height and the distance from the bound 
state to the barrier x*(t). The simplest stochastic prop- 
erties of them are listed below: 



(AG*(t)) = AG* 
<**(*)> =4 

(z*(t)xt(0)) - (^(0)) 2 = K x (t) 
(AG*(t)AG*(0)) - (AG*(0)) 2 = K g (t) 
(AG*(i)*t(0)) - (AG*(0))(^(0)) = K gx (t). 



(5) 



Here a stationary process and finite time correlation func- 
tions are assumed. Then using cumulant expansion of 
(kf(t)} again and truncating it to second order, we have 



(kf(t)} = fc exp 



-/3AG* + ^K g - M. 



f3K gx ) 2 



x exp 



(3 2 K X 



PK g 



2K : 
T 



(3K X 



(6) 



where K g , K x , and K gx are the variance and covariance 
of and between the two stochastic variables at the same 
time point. Because the above consideration has a simi- 
lar spirit with the Kubo- Anderson's stochastic line-shape 
theory [19| , and the second order truncations in variables 



Different from the Bell expression, when the force is 
larger, the dissociate rate increases exponentially with 
the square of force. This conclusion is very similar with 
that of Dembo et al. [^| . However the physical origin is 
completely different: the square of the force here arises 
from the fluctuation of the distance x* . It also means the 
bond is still slip and the lifetime of the bond is shorter 
than that predicted by the Bell formula. Of course, if 
the variance of the distance K x is very small, the mod- 
ification to Bell expression can be neglected. Although 
the result is interesting, in the following part we only fo- 
cus on the fourth case, in which (iv) both the distance 
and the barrier height are stochastic variables. Because 
the force is positive at the beginning, if K gx < or 
(l3K gx — xl) < 0, the behavior of the average dissoci- 
ation rate is then similar to the case in (iii). However if 



rt = 



/3K gx 



> 0, we see that the rate first decreases 



with the increasing of the force, and then increases when 
the force is beyond f c = x\/ j3K x , where the new parame- 
ters x\ and f c are defined for they have same dimensions 
of Distance and Force. Hence our model predicts the pos- 
sibility of the catch-slip bond transition at some critical 
transition force f c . 

We first consider the single molecule constant force 
rupture experiment [Toj |. where the average lifetime of 
the bond sPSGL-1— P-selectin was measured. Because 
the average dissociation rate kf(t) is time independent, 
the survival probability P(t) is a simple exponential func- 
tion exp[— t(kf)]. Following the general definition, the 
average lifetime of the bond is just t(f) = 1/ (kf). There 
are six parameters in this model. But in fact we can 
combine them into only three: the intrinsic dissociation 



rate k^ 



ko exp 



-0AGl 
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which is the param- 



eter that experiments can measure in practice, and the 
effective distance x\ defined above, and K x ; they are in- 
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dependent of each other. The average time then is 




(8) 



It is unexpected to find that the average lifetime of the 
bond is a Gaussian-like function with respect to the force: 
the mean value is f c , the variance a and a prefactor 
TV -1 ; their corresponding definitions see the above equa- 
tion. Apparently they are still independent. Accord- 
ing to the characteristic of Gaussian function, we can 
easily estimate the relevant parameters from the experi- 
mental data even without numerical methods; see Fig.^ 
they are respectively Hq « 133.0 /s, x\ « 2.88 nm, and 
K x w 1 nm 2 , whereas the important correlation coeffi- 
cient K gx > 12.0 pN • nm 2 and the catch-slip transition 
force / c w 12 pN which is directly read out from the 
experimental data. Here the estimation is performed at 
room temperature. We see that our prediction agrees 
with the data very well, in particular when the force is 
lower than f c 20]. Interestingly, the above result also 
shows that, only through forced dissociation experiment, 
we cannot isolate the precise information about the vari- 
ance of the energy barrier height and the cross variance 
of the height and the distance. 




OPSGL-l 
□ sPSGL-1 
— GSRM by us 



□ 

□ 



20 30 40 

Force (pN) 

FIG. 1: (Color online) Average lifetime as a function of 
the applied constant force for bonds of dimeric P-selectin 
with monomeric sPSGL-1 (blue square symbols) |1C| and the 
rescaled dimeric PSGL-1 (red circle symbols) from Ref. (l4|. 
The blue solid line is given by GSRM. The definitions of the 
symbols N, u, and f c see Eq. JSJ. 



More challenging experiment to our theory is force 
jump-ramp case, where the force increases linearly in 



time from a initial jump force fo, / = fo + rt, and r 
is ramp rate 0]. -"- n general, the stationary assumption 
of the forced dissociation processes should be more rea- 
sonable for the constant case. Considering that the direct 
extension of the Bell expression to time-dependent force 
case still provides insightful results , it is of interest to 
see what we can get by extending Eq. (JSJ) to force jump- 
ramp cases. Because the experimental data is typically 
presented in terms of the force histogram, we calculate 
the rupture force distribution P(f, fo) according to defi- 
nition P(f, fo)df = —(dS/dt)dt, which is given below, 



P(/,/o) = -exp 
r 



if - fc) 2 N ft (f-M* 

— / df e *>* 

2(7 r Jfo 



(9) 



We see that the average lifetime can be extracted from 
the above equation by setting / = / , i.e., t(f) = 
l/rP(/o,/o). 

We calculate the force distributions of the steady 
ramps (fo = pN) at ramp rates 210 and 1400 pN/s 
to compare with the experiment performed by Evans et 
al. 0]; see Fig. El where we use the same parameters 
obtained above. We find that the main qualitative char- 
acteristics of the predictions and the experimental data 
are the same: the distributions and the force histograms 
reach the maximum and minimum at two distinct forces, 
which are named / m j n and / max respectively. This obser- 
vation could be understood by setting the derivative of 
Eq. ijjjj with respect to / equal to zero, 



, , Aa 2 
./ - fc = exp 



(f-fc) 2 
2<r 2 



(10) 



Interestingly, Eq. (|10fl has no solutions when the loading 
rate is smaller than a critical ramp rate r c , which can be 
obtained by simultaneously solving the above equation 
and its first derivative 



1 



N 



(/ - fc) exp 



(/ - fc) 2 



2ct 2 



(11) 



Then r c = N(f* ~ fc) exp[(/* - f c ) 2 /2a% here /* is the 
force at which the maximum and minimum merge. We 
estimate r c w 9 pN/s using the current parameters. If 
f < f c , then the distribution is monotonous and decreas- 
ing function. Another important prediction is that the 
values of / m i n and / max must be larger than the catch-slip 
transition force f c . Indeed the experimental observation 
shows that the force values at the minimum force his- 
togram are around a certain values even the ramp rates 
change 10 folds, (see Figs. 2 and 4 in Ref. H3)- Accord- 
ing to Eq. H10JI . when the ramp rate is sufficiently large, 
we easily obtain 



fn 



fc 



No- 
r 



(12) 



by linear expansion. Therefore we predict that / m i n s ob- 
served in Evans et al. experiment |l2j | are almost the 
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catch-slip transition force observed in the constant force 
rupture experiment performed by Marshall et al. [Ioll2l| . 
Unfortunately, simply analytic relationship between / max 
and r cannot be found from Eq. flS) . Even so, the 
extrema equation implies that / max is a monotonous 
and increasing function of the ramp r, but the increas- 
ing is very slow and is about / ma x oc y/lnr instead of 
/max oc In r 
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FIG. 2: (Color online) The distribution of rupture forces un- 
der two loading rates predicted by our model for binding of 
P-selectin to sPSGL-1. The symbols of inset are the steady 
ramp experimental data [12|, where the loading rates are 210 
pN/s (red squares) and 1400 pN/s (blue circles), respectively. 
We must point out that the experimental data is for binding of 
P-selectin to PSGL-1, while our parameters are for sPSGL-1. 



We know that the rupture force distribution of simple 
slip bond only has a maximum at a certain force value 
that depends on the ramp rate j^. Therefor the catch- 
slip bond can easily be distinguished from the slip case 
by the presence of a minimum on the force rupture den- 
sity function at nonvanishcd force. Because the above 
analysis is independent of the initial force jump /q, in or- 
der to track the catch behaviors in the force jump-ramp 
experiments, the initial force /o should be smaller than 

fc 

In conclusion, we proposed a stochastic dissociation 
rate model to explain the intriguing catch-slip bond tran- 
sitions observed in the single molecule forced dissociation 
experiments, while the fluctuating rate is dependent on 
the two correlated stochastic control variables, the en- 
ergy barrier height AG* and the distance x* between 
the bound state and the energy barrier. Compared to 
the previous models with five [l^ . seven and four 
parameters |l4j involved, our model only requires three 
physical parameters: fcp, x\, and K x . Moreover, GSRM 
does not need the classical pathway concept or a priori 



catch bond 0, 0, Q, |2^| . Because there is no direct 
experiments or molecular structures supporting the path 
way scenario, a change in concept would be important 
for further experimental study of the catch-slip bonds. 

F.L. would like to thank Tokyo Institute of Technology 
for hospitality, where this work was completed. 
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